A New Constitutive Model Based on Taylor Series and Partial Derivatives for Predicting High-Temperature Flow Behavior of a Nickel-Based Superalloy

Ni-based superalloys are widely used in aerospace applications. However, traditional constitutive equations often lack the necessary accuracy to predict their high-temperature behavior. A novel constitutive model, utilizing Taylor series expansions and partial derivatives, is proposed to predict the high-temperature flow behavior of a nickel-based superalloy. Hot compression tests were conducted at various strain rates (0.01 s−1, 0.1 s−1, 1 s−1, and 10 s−1) and temperatures (850 °C to 1200 °C) to gather comprehensive experimental data. The performance of the new model was evaluated against classical models, specifically the Arrhenius and Hensel–Spittel (HS) models, using metrics such as the correlation coefficient (R), root mean square error (RMSE), sum of squared errors (SSE), and sum of absolute errors (SAE). The key findings reveal that the new model achieves superior prediction accuracy with an R value of 0.9948 and significantly lower RMSE (22.5), SSE (16,356), and SAE (5561 MPa) compared to the Arrhenius and HS models. Additionally, the stability of the first-order partial derivative of logarithmic stress with respect to temperature (∂lnσ/∂T) indicates that the logarithmic stress–temperature relationship can be approximated by a linear function with minimal curvature, which is effectively described by a second-degree polynomial. Furthermore, the relationship between logarithmic stress and logarithmic strain rate (∂lnσ/∂lnε˙) is more precisely captured using a third-degree polynomial. The accuracy of the new model provides an analytical basis for finite element simulation software. This helps better control and optimize processes, thus improving manufacturing efficiency and product quality. This study enables the optimization of high-temperature forming processes for current superalloy products, especially in aerospace engineering and materials science. It also provides a reference for future research on constitutive models and high-temperature material behavior in various industrial applications.


Introduction
Superalloys, particularly those based on nickel, are widely used in aerospace, power generation, and various high-temperature settings due to their remarkable mechanical characteristics and ability to withstand harsh conditions [1][2][3].These alloys are designed to endure high temperatures, corrosive environments, and significant mechanical stresses, making them ideal for gas turbine engines, aircraft structures, and nuclear reactors [4,5].Among the nickel-based superalloys, GH4169, also known as Inconel 718, has gained considerable attention for its outstanding creep resistance, high strength, and exceptional corrosion resistance at high temperatures [6][7][8].The material used in this paper is a typical GH4169 superalloy.The GH4169 superalloy is widely used in modern aero engines, gas turbines, and high-temperature components [9].Its importance has led researchers to study the relationship between its plastic flow stress and deformation conditions.Constitutive models play a crucial role in understanding and predicting the flow behavior of materials under various processing conditions, such as temperature, strain rate, and strain [10].A precise constitutive model is essential for optimizing the manufacturing processes of GH4169 superalloy components, such as hot forging, extrusion, and rolling [11].By predicting flow stress under varying processing conditions, the model helps engineers design and control the manufacturing processes more effectively, resulting in enhanced product quality and lower production costs [12].
Classical constitutive models, such as the Arrhenius and Hensel-Spittel (HS) models, have been widely used to describe the high-temperature deformation behavior of various materials [13,14].However, these models often lack the accuracy to effectively simulate the complex deformation behaviors of GH4169 under varying operational conditions.The Arrhenius model, based on the theory of thermally activated processes, presumes a linear relationship between the logarithm of flow stress and the logarithm of strain rate [15].In contrast, the HS model is phenomenological, considering the influence of temperature, strain rate, and strain on flow stress through a multiplicative formula [16].Despite the extensive use, the accuracy of these classical models in predicting the high-temperature flow behavior of the GH4169 superalloy has been questioned [17].Numerous studies have highlighted discrepancies between predicted and experimental flow stress values, particularly at elevated strain rates and temperatures [18,19].Lin et al. [20] developed models considering the coupled effects of deformation temperatures and strain rate.Xiao et al. [21] established a new physic-based constitutive model for GH4169 based on experimental results, considering the effects of dislocation movement.Zhou et al. [22] developed constitutive models for GH4169 that predict flow stresses considering the coupled effects of deformation parameters during work hardening-dynamic recovery (DRV) and dynamic recrystallization (DRX) periods.
In recent years, advanced mathematical techniques like artificial neural networks (ANNs) and support vector regression (SVR) have been applied to develop constitutive models for various materials [23,24].These methods have shown promising results in capturing the complex nonlinear relationships between flow stress and deformation variables.For instance, Zhu et al. [25] developed an ANN-based constitutive model for the GH4169 superalloy, which demonstrated improved accuracy compared to classical models.Similarly, Wen et al. [26] employed SVR to predict the flow behavior of the GH4169 superalloy, achieving high prediction accuracy.Zheng et al. [27] established a constitutive model for GH4169 based on the Arrhenius equation and optimized the model parameters using a backpropagation artificial neural network (BP-ANN).Gu et al. [28] constructed a data-driven constitutive model based on experimental data to describe the deformation behavior of the GH4169 superalloy.Wen et al. [29] proposed a two-stage deep learningbased constitutive model framework for GH4169 that incorporates a deformation-induced temperature rise and its influence on flow stress, comparing it favorably to the traditional Arrhenius-type model.However, the interpretability of these models remains a challenge, and their application in the GH4169 superalloy is somewhat limited.The models presented in these papers are difficult to apply to numerical simulation software, thus highlighting the need to develop a new model with a mathematical expression that can be used in numerical simulation software.
Recent advancements have highlighted the potential of higher-order mathematical approaches to improve the accuracy of constitutive models.The use of Taylor series expansions and partial derivatives in constitutive modeling has attracted considerable attention.These approaches offer a promising path for developing more accurate predictive models by capturing subtle variations in material behavior under different conditions [30].In the context of the GH4169 superalloy, limited research has been conducted on the application of Taylor series expansions and partial derivatives for constitutive modeling.This study aims to address these gaps by proposing a novel constitutive model that leverages Taylor series expansions and partial derivatives to predict the high-temperature flow behavior of GH4169 superalloy.By conducting comprehensive hot compression tests and evaluating the new model's performance against classical models, this study seeks to provide a more accurate and reliable tool for predicting the flow stress of GH4169 under various strain rates and temperatures.The development of this model will advance constitutive modeling techniques and optimize manufacturing processes for GH4169 superalloy components in aerospace and other high-temperature industries.

Material
The experimental material used in this study is the wrought GH4169 superalloy, supplied by Deyang Wanhang Die Forging Co., Ltd. in Deyang, China, and its nominal chemical composition is listed in Table 1.The alloying elements in GH4169 are Ni, Fe, Cr, Nb, and Mo.The original microstructure of the wrought GH4169 alloy is shown in Figure 1.attention.These approaches offer a promising path for developing more accurate predic tive models by capturing subtle variations in material behavior under different conditions [30].In the context of the GH4169 superalloy, limited research has been conducted on the application of Taylor series expansions and partial derivatives for constitutive modeling This study aims to address these gaps by proposing a novel constitutive model that lever ages Taylor series expansions and partial derivatives to predict the high-temperature flow behavior of GH4169 superalloy.By conducting comprehensive hot compression tests and evaluating the new model's performance against classical models, this study seeks to pro vide a more accurate and reliable tool for predicting the flow stress of GH4169 under var ious strain rates and temperatures.The development of this model will advance constitu tive modeling techniques and optimize manufacturing processes for GH4169 superalloy components in aerospace and other high-temperature industries.

Material
The experimental material used in this study is the wrought GH4169 superalloy, sup plied by Deyang Wanhang Die Forging Co., Ltd. in Deyang, China, and its nominal chem ical composition is listed in Table 1.The alloying elements in GH4169 are Ni, Fe, Cr, Nb and Mo.The original microstructure of the wrought GH4169 alloy is shown in Figure 1.

Experimental Procedures and Results
Hot forming is a critical technology for superalloys in the aerospace industry.Ele vated temperatures during hot forming enable significant plastic deformation of metallic alloys that exhibit limited ductility at room temperature.Hot forging processes for super alloys must achieve optimal high-temperature mechanical properties to meet perfor mance specifications under extreme thermal conditions [31].To investigate the hot form ing of the GH4169 superalloy in detail, eight temperatures (850 °C, 900 °C, 950 °C, 1000 °C, 1050 °C, 1100 °C, 1150 °C, and 1200 °C) and four strain rates (0.01 −1 , 0.1 −1 , 1 −1 , and 10 −1 were selected.Consequently, 32 cylindrical samples were processed from the wrough rod.Hot deformation samples with a height of 9 mm and a diameter of 6 mm were pre pared using wire electrodischarge machining.The samples were individually heated to the target temperature inside the Gleeble-3500 thermal compression testing machine (Data Sciences International, Inc., St. Paul, MN, USA) at a rate of 10 °C/s.Subsequently they were isothermally compressed at different strain rates (0.01 −1 , 0.1 −1 , 1 −1 , and 10 −1 ) to achieve a deformation of 60%.The isothermal compression process and corresponding

Experimental Procedures and Results
Hot forming is a critical technology for superalloys in the aerospace industry.Elevated temperatures during hot forming enable significant plastic deformation of metallic alloys that exhibit limited ductility at room temperature.Hot forging processes for superalloys must achieve optimal high-temperature mechanical properties to meet performance specifications under extreme thermal conditions [31].To investigate the hot forming of the GH4169 superalloy in detail, eight temperatures (850 C, and 1200 • C) and four strain rates (0.01 −1 , 0.1 −1 , 1 −1 , and 10 −1 ) were selected.Consequently, 32 cylindrical samples were processed from the wrought rod.Hot deformation samples with a height of 9 mm and a diameter of 6 mm were prepared using wire electrodischarge machining.The samples were individually heated to the target temperature inside the Gleeble-3500 thermal compression testing machine (Data Sciences International, Inc., St. Paul, MN, USA) at a rate of 10 • C/s.Subsequently, they were isothermally compressed at different strain rates (0.01 −1 , 0.1 −1 , 1 −1 , and 10 −1 ) to achieve a deformation of 60%.The isothermal compression process and corresponding experimental parameters are shown in Figure 2 and Table 2, respectively.After each sample completed the isothermal compression test, it was quenched in water.The flow stress-strain curves were automatically recorded from the load-displacement data during the isothermal compression experiments.experimental parameters are shown in Figure 2 and Table 2, respectively.After each sample completed the isothermal compression test, it was quenched in water.The flow stressstrain curves were automatically recorded from the load-displacement data during the isothermal compression experiments.Figure 3 shows the stress-strain curves at different strain rates and deformation temperatures.At the same strain rate, the stress gradually decreases with increasing deformation temperature, and the peak stress also decreases as the temperature rises.For a given strain rate at a specific temperature, the stress sharply increases with strain during the initial stages of deformation until it reaches its maximum value.This increase is primarily attributed to work hardening.As deformation continues, dislocation slip occurs, leading to a sharp increase in dislocation density within the grains, resulting in dislocation interactions.In the later stages, stress decreases with increasing strain due to dynamic recovery and recrystallization.When the deformation reaches a certain level, the stored energy becomes the driving force for recrystallization, initiating dynamic recrystallization.The occurrence of dynamic recrystallization leads to the disappearance of internal dislocations, reducing the material's deformation resistance.As deformation progresses, the recrystallization continues until it is complete.After complete internal dynamic recrystallization, the flow stress reaches a steady-state stage.Furthermore, the non-monotonic behavior in the plot at 1200 °C in Figure 3d (depicted by the blue curve at the bottom) shows an initial increase in stress with strain, reaching a peak stress at a relatively low strain value (around 0.05).After reaching this peak, the stress decreases significantly, exhibiting a softening behavior.The stress continues to drop as the strain increases, indicating dynamic recrystallization or flow softening mechanisms at this elevated temperature.This behavior suggests that at 1200 °C, the material undergoes significant softening after the initial hardening, likely due to the thermal activation of mechanisms that reduce  Figure 3 shows the stress-strain curves at different strain rates and deformation temperatures.At the same strain rate, the stress gradually decreases with increasing deformation temperature, and the peak stress also decreases as the temperature rises.For a given strain rate at a specific temperature, the stress sharply increases with strain during the initial stages of deformation until it reaches its maximum value.This increase is primarily attributed to work hardening.As deformation continues, dislocation slip occurs, leading to a sharp increase in dislocation density within the grains, resulting in dislocation interactions.In the later stages, stress decreases with increasing strain due to dynamic recovery and recrystallization.When the deformation reaches a certain level, the stored energy becomes the driving force for recrystallization, initiating dynamic recrystallization.The occurrence of dynamic recrystallization leads to the disappearance of internal dislocations, reducing the material's deformation resistance.As deformation progresses, the recrystallization continues until it is complete.After complete internal dynamic recrystallization, the flow stress reaches a steady-state stage.Furthermore, the non-monotonic behavior in the plot at 1200 • C in Figure 3d (depicted by the blue curve at the bottom) shows an initial increase in stress with strain, reaching a peak stress at a relatively low strain value (around 0.05).After reaching this peak, the stress decreases significantly, exhibiting a softening behavior.The stress continues to drop as the strain increases, indicating dynamic recrystallization or flow softening mechanisms at this elevated temperature.This behavior suggests that at 1200 • C, the material undergoes significant softening after the initial hardening, likely due to the thermal activation of mechanisms that reduce dislocation density or promote grain boundary sliding.To investigate the constitutive behavior of the superalloy, the flow stress curve was discretized into ten equal intervals, covering a strain range from 0.04 to 0.9.The stress values in Table 3 were obtained by linear interpolation of the flow curves shown in Figure 3.The colors in Table 3 transition from light to dark, from white to red, where a deeper red background indicates higher stress values.Analyzing the data for a specific strain shows that, at the same temperature, the stress value increases with the strain rate.Conversely, at the same strain rate, the stress value gradually decreases with an increase in temperature.Comparing these 10 sets of data clearly shows that they exhibit the same trend.from light to dark, from white to red, where a deeper red background indicates high stress values.Analyzing the data for a specific strain shows that, at the same temperatu the stress value increases with the strain rate.Conversely, at the same strain rate, the stre value gradually decreases with an increase in temperature.Comparing these 10 sets data clearly shows that they exhibit the same trend.Table 3. Stress matrix corresponding to different strain levels (the units of strain rate, temperatu and stress are s −1 , K, and MPa, respectively).

Classical Constitutive Model
For accurate numerical simulations, a reliable flow curve prediction model must be established.The flow curve describes the relationship between the flow stress and strain of a material at different temperatures and strain rates.A constitutive model is a mathematical expression used to fit the experimental data of the flow curve.In this paper, a comparison is made between two commonly used constitutive models, the Arrhenius model and the Hensel-Spittel (HS) model, before proposing a new phenomenological constitutive model.The Arrhenius model is based on activation energy theory, which assumes that the logarithm of flow stress has a linear relationship with temperature and strain rate.The Hensel-Spittel (HS) model is an empirical model that assumes flow stress has a linear relationship with a power function of temperature and strain rate.

Arrhenius Model
The Arrhenius model was first proposed by Sellars and Tegart [32].This widely used constitutive model describes the relationship between flow stress, temperature, and strain rate during high-temperature plastic deformation of alloys.It assumes that flow stress is a hyperbolic sine function of temperature T and deformation activation energy Q, with its mathematical expression outlined as follows: .
where .ε is the strain rate of the material (s −1 ), σ is flow stress of the material (MPa), Q is the activation energy of thermal deformation (J•mol −1 ), R is the universal gas constant (8.314J•K −1 •mol −1 ), T is the thermodynamic temperature of the material (K), α is the stress level parameter, A is the structure factor (s −1 ), n is the stress index, and A, α, n, β (β = αn) and Q are the material parameters.
Taking the logarithm of the Equation (1): Multivariate nonlinear regression based on Equation ( 2) was performed with the flow stress data to solve for the material parameters lnA, α, n and Q at different strain levels (Table 3).Then, the curves of lnA-ε, α-ε, n-ε, and Q-ε (as shown in Figure 4) were fitted with polynomials to obtain the constitutive equation with strain compensation.
Materials 2024, 17, x FOR PEER REVIEW 7 of 25 (1) where is the strain rate of the material (s −1 ),  is flow stress of the material (MPa),  is the activation energy of thermal deformation (J•mol −1 ),  is the universal gas constant (8.314J•K −1 •mol −1 ),  is the thermodynamic temperature of the material (K),  is the stress level parameter,  is the structure factor (s −1 ),  is the stress index, and , , ,  ( = ) and  are the material parameters.
Taking the logarithm of the Equation ( 1): (2) Multivariate nonlinear regression based on Equation ( 2) was performed with the flow stress data to solve for the material parameters ln, ,  and  at different strain levels (Table 3).Then, the curves of ln-, -, -, and - (as shown in Figure 4) were fitted with polynomials to obtain the constitutive equation with strain compensation.As shown in Figure 4a, the regression values of ln start at around 40 and rise sharply to around 55 at a strain level of 0.2.After reaching this peak, the values decrease slightly and then stabilize around 55 as the strain increases further, exhibiting a gentle upward trend toward a strain level of 0.9.As shown in Figure 4b, the regression values of  start at 0.005 and decrease rapidly to about 0.003 at a strain level of 0.2.After this initial As shown in Figure 4a, the regression values of lnA start at around 40 and rise sharply to around 55 at a strain level of 0.2.After reaching this peak, the values decrease slightly and then stabilize around 55 as the strain increases further, exhibiting a gentle upward trend toward a strain level of 0.9.As shown in Figure 4b, the regression values of α start at 0.005 and decrease rapidly to about 0.003 at a strain level of 0.2.After this initial drop, the values gradually increase, showing a slight upward trend toward a strain level of 0.9.As shown in Figure 4c, the regression values of n start at about 6 and decrease significantly to about 4.5 at a strain level of 0.2.After this initial drop, the values stabilize around 4.5, showing a gentle downward trend as the strain increases further.As shown in Figure 4d, the regression values of Q start at about 5.5 and rise sharply to about 7 at a strain level of 0.2.After reaching this peak, the values decrease slightly and then stabilize around 7 as the strain increases further, exhibiting a gentle upward trend toward a strain level of 0.9.
A sixth-order polynomial is used to fit these data, thereby obtaining parametercontaining expressions as shown in Equation (3).The polynomial order is determined by the regression accuracy, but a higher order may cause overfitting.
By substituting Equation (3) into Equation ( 1), the Arrhenius constitutive equation for GH4169 alloy can be obtained.To visually observe the prediction accuracy of the model, the experimental data and predicted data from the Arrhenius constitutive equation are plotted on the same graph for comparison, as shown in Figure 5.
drop, the values gradually increase, showing a slight upward trend toward a strain l of 0.9.As shown in Figure 4c, the regression values of n start at about 6 and decr significantly to about 4.5 at a strain level of 0.2.After this initial drop, the values stab around 4.5, showing a gentle downward trend as the strain increases further.As sh in Figure 4d, the regression values of  start at about 5.5 and rise sharply to about 7 strain level of 0.2.After reaching this peak, the values decrease slightly and then stab around 7 as the strain increases further, exhibiting a gentle upward trend toward a st level of 0.9.
A sixth-order polynomial is used to fit these data, thereby obtaining parametertaining expressions as shown in Equation (3).The polynomial order is determined by regression accuracy, but a higher order may cause overfitting.By substituting Equation (3) into Equation ( 1), the Arrhenius constitutive equa for GH4169 alloy can be obtained.To visually observe the prediction accuracy of model, the experimental data and predicted data from the Arrhenius constitutive equa are plotted on the same graph for comparison, as shown in Figure 5. Comparing the curves at four different strain rates reveals that as the strain rat creases, the material's flow stress level increases significantly.This indicates the strain sensitivity of the material.At high strain rates, dislocation movement is impeded, re ing in significant material strengthening.At a given strain rate, as the temperatur creases, the flow stress gradually decreases.At high temperatures, thermal activation Comparing the curves at four different strain rates reveals that as the strain rate increases, the material's flow stress level increases significantly.This indicates the strain rate sensitivity of the material.At high strain rates, dislocation movement is impeded, resulting in significant material strengthening.At a given strain rate, as the temperature increases, the flow stress gradually decreases.At high temperatures, thermal activation processes intensify, dislocation motion resistance decreases, and the material softens.The flow stress peaks at 850 • C and is lowest at 1200 • C. In the low strain range, the model's predicted values slightly deviate from the experimental values.It is speculated that in the small strain stage, complex microstructural evolution behaviors such as strain hardening and dynamic recovery may occur, which the Arrhenius model cannot fully describe.At a high strain rate of 10 s −1 , the model's deviation slightly increases.This suggests that the Arrhenius model is more suitable for low strain rate conditions.At high strain rates, the introduction of other dynamic softening factors may be necessary.Overall, the Arrhenius model can describe high-temperature rheology to some extent, but its quantitative accuracy requires further improvement.

Hensel-Spittel Constitutive Model
The Hensel-Spittel (HS) constitutive model accurately describes the relationship between flow stress, strain, strain rate, and temperature during the thermal deformation of metallic materials [33].By selecting appropriate model parameters, the HS constitutive model can predict the mechanical behavior of various metallic materials under different processing parameters, providing an important basis for optimizing and controlling the hot forming process [34].Furthermore, the HS constitutive model has gained extensive application in hot forming numerical simulations due to its simple mathematical form and easily accessible material parameters [35][36][37][38].Renowned commercial finite element software like Forge NxT has incorporated this model into their material constitutive libraries, showcasing its universality and effectiveness in the engineering field, as shown in Equation ( 4): where ε, σ, .
ε, and T represent strain, flow stress, strain rate, and temperature, respectively, and A and m 1 ~m8 are material constants.Taking the natural logarithm of Equation ( 4): There is a linear relationship between lnσ and T, lnε, ln .
ε, 1/ε, Tln(1 + ε), ε, Tln .ε, and Tlnε.Therefore, solving for the material constants in Equation ( 5) is a typical multivariate linear regression problem.Similarly, the strain is divided into ten equal intervals between 0.04 and 0.9.Then, by interpolating the original compression data (as shown in Figure 3), the flow stress at each strain level can be obtained.Using MATLAB software (Version R2016b) to solve for the parameters, the multivariate linear regression results for the material parameters are listed in Table 4.The comparison of the HS model's predicted flow stress with the experimental flow stress is shown in Figure 6.

Evaluation of the Applicability of Traditional Constitutive Models
To visually evaluate the predictive capability of the Arrhenius model and the Hensel Spittel model for the high-temperature flow behavior of wrought GH4169 alloy, their cal culated values were compared with experimental values.Figures 5 and 6, respectively show the comparison of the predicted and experimental flow stress values for the Arrhe nius model and the HS model.It can be seen that the error is larger at lower temperature of 850 °C, 900 °C, and 950 °C.As shown in Figure 6, the predictive accuracy of the H constitutive equation is lower at lower temperatures (850 °C, 900 °C, and 950 °C), while i is higher in the high-temperature range (1000 °C, 1050 °C, 1100 °C, 1150 °C, and 1200 °C) In conclusion, the traditional Arrhenius and Hensel-Spittel models can describe th high-temperature flow behavior of the wrought GH4169 superalloy to some extent.Over all, the predictive accuracy of the HS model is higher than that of the Arrhenius mode but both models still face challenges in further improving their predictive accuracy.Th main reasons for the prediction deviations are that neither model considers the depend ence of flow stress on strain, whereas the flow stress of the GH4169 superalloy actuall changes significantly with strain.Additionally, traditional constitutive models lack con sideration of the interaction between different deformation parameters (such as tempera ture, strain rate, and strain).Therefore, it is necessary to develop new constitutive model to more accurately describe the high-temperature flow characteristics of the wrough GH4169 superalloy.

Evaluation of the Applicability of Traditional Constitutive Models
To visually evaluate the predictive capability of the Arrhenius model and the Hensel-Spittel model for the high-temperature flow behavior of wrought GH4169 alloy, their calculated values were compared with experimental values.Figures 5 and 6 In conclusion, the traditional Arrhenius and Hensel-Spittel models can describe the high-temperature flow behavior of the wrought GH4169 superalloy to some extent.Overall, the predictive accuracy of the HS model is higher than that of the Arrhenius model, but both models still face challenges in further improving their predictive accuracy.The main reasons for the prediction deviations are that neither model considers the dependence of flow stress on strain, whereas the flow stress of the GH4169 superalloy actually changes significantly with strain.Additionally, traditional constitutive models lack consideration of the interaction between different deformation parameters (such as temperature, strain rate, and strain).Therefore, it is necessary to develop new constitutive models to more accurately describe the high-temperature flow characteristics of the wrought GH4169 superalloy.

Development of a New Constitutive Model
The Arrhenius model highlights that the key to constructing the constitutive equation under specific strain is determining the mathematical relationship among stress, strain rate, and temperature.However, the Arrhenius model itself is an implicit equation, making it relatively complex to solve in engineering applications.Additionally, the fitting accuracy of this model for the stress-strain values of the GH4169 alloy requires improvement.Theoretically, the relationship between stress, strain rate, and temperature can be represented as σ = f ε, .ε, T .However, due to the significant nonlinearity among these three variables, this study did not directly use this functional form.Before constructing a new constitutive model, the influence of strain rate and temperature on stress was investigated.To this end, experimental data at low strain (0.136), medium strain (0.518), and high strain (0.9) levels from Table 3 were analyzed.To reveal the intrinsic relationship between logarithmic stress, logarithmic strain rate, and temperature, the partial derivatives of logarithmic stress with respect to temperature and logarithmic strain rate were calculated using numerical discretization methods, such as finite differences.Common numerical discretization methods include forward difference, backward difference, and central difference.To ensure computational accuracy and stability, forward and backward differences were used at boundary points, while central differences were used at internal nodes.

Mathematical Principles
Partial derivatives are an important concept in multivariable calculus, indicating the rate of change in a function in the direction of a specific variable.For a two-variable function f (x, y), the partial derivative in the x-direction is defined as follows [39]: Similarly, the partial derivative of the function in the y-direction is defined as follows: In practical calculations of partial derivatives, numerical discretization methods are commonly employed, approximating derivatives using finite differences.Similarly, secondorder and higher-order derivatives can be obtained through the recursive application of first-order finite differences.To ensure accuracy when calculating derivatives at boundary points, one-sided difference schemes are typically used, whereas central difference schemes are often employed at internal nodes to balance accuracy and stability.
The Taylor series is an important local expansion of functions used to approximate complex nonlinear functions.If a function f (x) has continuous derivatives up to the nth order at the point x 0 , its Taylor expansion is given as follows [40]: Nickel-based superalloys exhibit significant nonlinear mechanical behavior under high strain rates and high-temperature conditions, with complex nonlinear coupling effects between stress levels, strain rates, and temperatures.Experimental data at strain levels of 0.136, 0.518, and 0.9 from Table 3 were selected for analysis.To quantitatively describe the functional relationship between logarithmic stress, temperature, and logarithmic strain rate, the partial derivatives of logarithmic stress with respect to temperature and logarithmic strain rate were calculated using discrete formulas.The discrete formulas include forward difference, backward difference, and central difference.Forward difference and backward difference were used at the initial and final boundaries, respectively, while central difference was used in the intermediate region.Based on the discrete difference method, the nth-order partial derivatives of logarithmic stress with respect to temperature and logarithmic strain rate at any point on the experimental stress-strain curve can be calculated using Formula (9) as the central difference.
where n represents the order of the partial derivative, and i and j represent the indices of strain rate and temperature, respectively.For example, when i = 1 and j = 1, lnσ ij denotes the logarithmic stress at a strain rate of 0.01 s −1 and a temperature of 850 • C. The maximum values of i and j correspond to the number of strain rates and temperatures, respectively.When i = 1 and j = 1, Formula (9) transforms into the backward difference Formula (10).
From Formulas ( 9)- (11), it is evident that the accuracy of calculating the partial derivatives of stress is closely related to the step size of the finite differences.The denser the experimental points for temperature and strain rate, the closer the numerical calculation results are to the actual partial derivatives.Maximizing the number of experimental data points is crucial for accurately describing the constitutive behavior of the wrought GH4169 superalloy.The first, second, and third-order partial derivatives of logarithmic stress with respect to temperature and logarithmic strain rate were calculated at strain levels of 0.136, 0.518, and 0.9.The relevant results are shown in Figures 7-9, respectively.For the three strain levels, logarithmic stress decreases with increasing deformation temperature (as shown in Figures 7a, 8a and 9a).This indicates that, at a given strain level, higher strain rates lead to higher stress values.At a specific strain level, the first-order partial derivatives of logarithmic stress with respect to temperature (∂lnσ/∂T) show little difference at different temperatures, with a magnitude of about 0.001, as shown in Figures 7b, 8b and 9b.This indicates that the rate of change in logarithmic stress with temperature (i.e., the slope of the logarithmic stress-temperature curve) remains relatively stable within the studied temperature range.The second-order partial derivatives (∂ 2 lnσ/∂T 2 ) and third-order partial derivatives (∂ 3 lnσ/∂T 3 ) of logarithmic stress with respect to temperature are close to zero at all strain levels, with smaller magnitudes, as shown in Figures 7c,d, 8c,d and 9c,d.
From the Taylor series Formula (8), it can be seen that if the Taylor series expansion of f (x) near the point x 0 is written out and the fourth and higher-order terms are ignored, the following can be obtained: The constant term f (x 0 ), which represents the value of the function at x 0 , determines the vertical position of the Taylor series expansion but does not affect the shape of the function.The first-order term is (x − x 0 ) f ′ (x 0 ), which is a linear function with a slope given by the first derivative f ′ (x 0 ).When x is close to x 0 , this term causes the function to exhibit approximately linear behavior near x 0 .The second-order term (x − x 0 ) 2 f ′′ (x 0 )/2! is a quadratic function with a coefficient determined by the second derivative f ′′ (x 0 ).If f ′′ (x 0 ) is not zero, this term introduces curvature, causing the function f (x) to deviate from linearity.However, if f ′′ (x 0 ) is close to zero, the contribution of this term becomes very small, and the function f(x) still exhibits approximately linear behavior near x 0 .The third-order term (x − x 0 ) 3 f ′′′ (x 0 )/3! is a cubic function with a coefficient determined by the third derivative f ′′′ (x 0 ).Similar to the second-order term, if f ′′′ (x 0 ) is close to zero, the contribution of this term also becomes very small.This implies that even if the function f (x) exhibits slight nonlinearity near x 0 , its rate of curvature change is very small, and it can still be well approximated by a linear function overall.
The first-order partial derivative of logarithmic stress with respect to temperature (∂lnσ/∂T) is relatively stable within the studied temperature range.This indicates that the logarithmic stress-temperature curve can be well approximated by a straight line with a constant slope within this temperature range.The second-order partial derivative (∂ 2 lnσ/∂T 2 ) and third-order partial derivative (∂ 3 lnσ/∂T 3 ) of logarithmic stress with respect to temperature are very close to zero within the studied temperature range.This suggests that the logarithmic stress-temperature curve has almost no curvature or rate of curvature change within this temperature range, exhibiting characteristics very close to linearity.However, to fit the experimental curve more precisely, higher-order partial derivatives are considered.Therefore, the relationship between logarithmic stress lnσ and temperature T within the studied temperature range can be described with high accuracy by a second-order function.
x FOR PEER REVIEW 13 of 25     From the Taylor series Formula ( 8), it can be seen that if the Taylor series expansion of () near the point  is written out and the fourth and higher-order terms are ignored, the following can be obtained: The importance of this conclusion is that it provides a simple and effective method for describing the law of variation of material flow stress with temperature.By analyzing the higher-order partial derivatives of logarithmic stress, the linearity of the lnσ-T relationship can be directly judged, thereby selecting an appropriate mathematical model to fit the experimental data.This not only greatly simplifies the process of establishing the constitutive model but also improves the computational efficiency and predictive capability of the model.It provides a theoretical foundation for establishing efficient and practical material constitutive models.This method is not only applicable to the high-temperature alloy studied in this paper but can also be extended to other types of materials, offering broad application prospects.
For the studied strain levels (low strain ε = 0.136, medium strain ε = 0.518, and high strain ε = 0.9), logarithmic stress increases with the increase in logarithmic strain rate (as shown in Figures 10a, 11a and 12a).The first-order partial derivative of logarithmic stress with respect to logarithmic strain rate (∂lnσ/∂ln .ε) varies significantly at all strain levels (as shown in Figures 10b, 11b and 12b), and the absolute value of the first-order partial derivative increases with the increase in temperature and logarithmic strain rate.This indicates that the impact of strain rate on flow stress becomes more significant with an increase in temperature and strain rate.In other words, the material is more sensitive to strain rate under high temperature and high strain rate conditions.Furthermore, it can be found that the second-order partial derivative of logarithmic stress with respect to logarithmic strain rate (∂ 2 lnσ/∂(ln .ε) 2 ) is relatively small at all strain levels (as shown in Figures 10c, 11c and 12c), while the third-order partial derivative (∂ 3 lnσ/∂(ln .ε) 3 ) is also relatively small (as shown in Figures 10d, 11d and 12d).To make the predicted values more accurately describe the experimental values, the relatively small third-order partial derivatives are also taken into account.Therefore, the relationship between logarithmic stress and logarithmic strain rate can be described by a cubic function.
temperature  within the studied temperature range can be described with high accuracy by a second-order function.
The importance of this conclusion is that it provides a simple and effective method for describing the law of variation of material flow stress with temperature.By analyzing the higher-order partial derivatives of logarithmic stress, the linearity of the ln- relationship can be directly judged, thereby selecting an appropriate mathematical model to fit the experimental data.This not only greatly simplifies the process of establishing the constitutive model but also improves the computational efficiency and predictive capability of the model.It provides a theoretical foundation for establishing efficient and practical material constitutive models.This method is not only applicable to the high-temperature alloy studied in this paper but can also be extended to other types of materials, offering broad application prospects.
For the studied strain levels (low strain  = 0.136, medium strain  = 0.518, and high strain  = 0.9), logarithmic stress increases with the increase in logarithmic strain rate (as shown in Figures 10a, 11a, and 12a).The first-order partial derivative of logarithmic stress with respect to logarithmic strain rate (ln ln ⁄ ) varies significantly at all strain levels (as shown in Figures 10b, 11b, and 12b), and the absolute value of the first-order partial derivative increases with the increase in temperature and logarithmic strain rate.This indicates that the impact of strain rate on flow stress becomes more significant with an increase in temperature and strain rate.In other words, the material is more sensitive to strain rate under high temperature and high strain rate conditions.Furthermore, it can be found that the second-order partial derivative of logarithmic stress with respect to logarithmic strain rate ( ln (ln ) ⁄ ) is relatively small at all strain levels (as shown in Fig- ures 10c, 11c, and 12c), while the third-order partial derivative ( ln (ln ) ⁄ ) is also relatively small (as shown in Figures 10d, 11d, and 12d).To make the predicted values more accurately describe the experimental values, the relatively small third-order partial derivatives are also taken into account.Therefore, the relationship between logarithmic stress and logarithmic strain rate can be described by a cubic function.

Establishment of the New Constitutive Model
Based on the previous analysis, the deformation of the constitutive equation with two

Establishment of the New Constitutive Model
Based on the previous analysis, the deformation of the constitutive equation with two

Establishment of the New Constitutive Model
Based on the previous analysis, the deformation of the constitutive equation with two variables can be expressed by the following formula: Then, using the Taylor series expansion and more accurately considering the coupling effect of temperature and strain rate, a new constitutive model is proposed, as shown in the following Formula ( 14): Further simplifying Formula (15): where ε represents the first-order response of stress to strain rate, describing the linear relationship between stress and logarithmic strain rate.The first-order term k 3 T represents the first-order response of logarithmic stress to temperature, reflecting the thermal softening effect of the material.The second-order term k 4 ln .ε 2 represents the nonlinear change in logarithmic stress with respect to the logarithmic strain rate, reflecting the nonlinear behavior of the material under strain rate.The cross term k 5 Tln . ε represents the interaction between strain rate and temperature on stress, describing how stress is influenced by both strain rate and temperature, and reflecting the coupling effect of the temperature and strain rate.The second-order term k 6 T 2 represents the nonlinear change in logarithmic stress with respect to temperature, reflecting the nonlinear characteristics of the material's thermal softening behavior.The third-order term k 7 ln .ε 3 represents the thirdorder response of stress to strain rate, also describing the complex nonlinear relationship between stress and logarithmic strain rate.The term k 8 ln .ε 2 T represents the interaction between the second-order term of strain rate and temperature, describing how stress is influenced by the square of logarithmic strain rate and temperature.The term k 9 ln .εT 2 represents the second-order interaction between strain rate and temperature, describing how stress is influenced by logarithmic strain rate and the square of temperature.The term k 10 ln .ε 3 T represents the interaction between the third-order term of the logarithmic strain rate and the first-order term of temperature, capturing the high-order nonlinear relationship between logarithmic strain rate and temperature.The term k 11 ln .εT 3 represents the interaction between the first-order term of the logarithmic strain rate and the thirdorder term of temperature.The term k 12 ln .ε 2 T 2 represents the interaction between the second-order term of the logarithmic strain rate and the second-order term of temperature, reflecting the more complex interaction between these two variables.
Based on the data provided in Table 3, the material parameters in Formula ( 15) can be obtained using multiple linear regression.Additionally, since different strains have different material parameters, the relationship between material parameters and strain can be fitted using a polynomial.The degree of the polynomial can be determined based on the regression accuracy.In this study, a seventh-degree polynomial was used to fit these material parameters, and the relationship between each material parameter and different strains was obtained, thus incorporating stress, strain, strain rate, and temperature into the new constitutive equation.The coefficients of each term are shown in Table 5, and the polynomial fitting curves of the quadratic model coefficients at different strain levels are shown in Figure 13.For different materials, the number of coefficients and the order of partial derivatives must be determined based on the complexity of their mechanical behavior.The accuracy of the new constitutive model is also ensured by these factors.Additionally, the coefficients must be solved mathematically.By substituting the relationship between the parameters of the new constitutive model for the GH4169 alloy and strain from Table 5 into Formula (15), a complete new constitutive model is obtained.This equation considers the nonlinear relationship between material parameters and strain, allowing the model to more accurately describe the mechanical behavior of the material under different deformation conditions.To evaluate the predictive ability of the newly established constitutive model, the stress predicted by the new model was compared with the experimentally measured data, as shown in Figure 14.In the figure, the solid lines of different colors represent the data obtained from experiments, reflecting the true mechanical behavior of the material under different temperatures and strain rates.The scatter points represent the stress values predicted by the new model under the corresponding conditions.By visually comparing the degree of agreement, the predictive accuracy and reliability of the new model can be qualitatively assessed.By substituting the relationship between the parameters of the new constitutive model for the GH4169 alloy and strain from Table 5 into Formula (15), a complete new constitutive model is obtained.This equation considers the nonlinear relationship between material parameters and strain, allowing the model to more accurately describe the mechanical behavior of the material under different deformation conditions.To evaluate the predictive ability of the newly established constitutive model, the stress predicted by the new model was compared with the experimentally measured data, as shown in Figure 14.In the figure, the solid lines of different colors represent the data obtained from experiments, reflecting the true mechanical behavior of the material under different temperatures and strain rates.The scatter points represent the stress values predicted by the new model under the corresponding conditions.By visually comparing the degree of agreement, the predictive accuracy and reliability of the new model can be qualitatively assessed.

Evaluation of Predictive Performance
To comprehensively evaluate the overall predictive accuracy of the classical model and the new model, a series of standard statistical parameters were introduced, including the correlation coefficient ( ), root mean square error (RMSE), sum of squared errors

Evaluation of Predictive Performance
To comprehensively evaluate the overall predictive accuracy of the classical model and the new model, a series of standard statistical parameters were introduced, including the correlation coefficient (R), root mean square error (RMSE), sum of squared errors (SSE), and sum of absolute errors (SAE).These parameters quantify the deviation between the model predictions and experimental measurements from different perspectives, providing an objective basis for comparing the performance of different models.The correlation coefficient (R) measures the linear correlation between the model predictions and the experimental values, and its calculation formula is as follows: where  17) and ( 18), respectively: where N represents the total number of data points.The smaller the values of RMSE and SSE, the smaller the deviation between the model predictions and the experimental values, indicating higher overall predictive accuracy of the model.The sum of absolute errors (SAE) calculates the cumulative amount of absolute errors between the model predictions and the experimental values.Its calculation formula is shown in Formula ( 19): This indicator is similar to RMSE and SSE.The smaller the value of SAE, the smaller the absolute deviation between the model predictions and the experimental values, indicating higher overall predictive accuracy of the model.
By calculating the values of these four statistical parameters for both the classical model and the new model, and comparing them, the overall predictive performance of the two models can be quantitatively evaluated.Generally, if the new model achieves better results in these parameters (e.g., R closer to 1, RMSE, SSE, and SAE smaller), it indicates that incorporating the nonlinear relationship between material parameters and strain effectively improves the overall predictive accuracy of the constitutive model.The quantitative comparison of the predictive accuracy of different models is shown in Table 6.
According to the quantitative comparison of the predictive abilities of different models in Table 6, the new model performs the best in all indicators.Its correlation coefficient (R) is 0.9948; the root mean square error (RMSE) is 22.5; the sum of squared errors (SSE) is 16,356; and the sum of absolute errors (SAE) is 5561 MPa.These values indicate that the model has superior fitting and predictive capabilities.The Arrhenius model follows closely, with a correlation coefficient (R) of 0.9888, a root mean square error (RMSE) of 44.1020, a sum of squared errors (SSE) of 19,450, and a sum of absolute errors (SAE) of 8832 MPa.However, its overall performance is slightly inferior to the new model.The HS model shows the poorest performance in all indicators, with a correlation coefficient (R) of 0.9827, a root mean square error (RMSE) of 44.0915, a sum of squared errors (SSE) of 38,881, and a sum of absolute errors (SAE) of 10,877 MPa, indicating its relatively weaker predictive ability.Figure 15 shows the relationship between the predicted values of different models and the experimental data.The solid line represents the experimental data, hollow circles represent the predictions of the HS model, plus signs represent the predictions of the Arrhenius model, and solid dots represent the predictions of the new model.By analyzing Figure 15a-d, it can be seen that under different strain rates and temperatures, the new model has the best fitting effect, accurately predicting the stress-strain relationship of the material under various conditions.The Arrhenius model follows, fitting the experimental data well under most conditions but performing slightly worse than the new model under high temperatures and high strain rates.The HS model has relatively poor predictive performance, especially under high temperature and high strain rate conditions, where its predictions deviate significantly from the experimental data.Therefore, the new model shows a clear advantage in predicting material behavior.Figure 15 shows the relationship between the predicted values of different models and the experimental data.The solid line represents the experimental data, hollow circles represent the predictions of the HS model, plus signs represent the predictions of the Arrhenius model, and solid dots represent the predictions of the new model.By analyzing Figure 15a-d, it can be seen that under different strain rates and temperatures, the new model has the best fitting effect, accurately predicting the stress-strain relationship of the material under various conditions.The Arrhenius model follows, fitting the experimental data well under most conditions but performing slightly worse than the new model under high temperatures and high strain rates.The HS model has relatively poor predictive performance, especially under high temperature and high strain rate conditions, where its predictions deviate significantly from the experimental data.Therefore, the new model shows a clear advantage in predicting material behavior.Figure 16 shows the predictive ability of different models for flow stress.The coefficient of determination R 2 for the Arrhenius model is 0.977; for the HS model, it is 0.967; and for the new model, it is 0.99.Combined with the previous analysis results, the new model has the best fitting effect under various strain rate and temperature conditions, accurately predicting the stress-strain relationship of the material and having the highest predictive accuracy.The Arrhenius model follows, with good predictive ability, but it is slightly inferior to the new model.The HS model has the poorest predictive ability, with significant deviations, especially under high-stress conditions.Therefore, the new model demonstrates the highest reliability and accuracy in practical applications.Consequently, this high-precision model, with its mathematical expression, can be integrated into forming simulation software.This integration provides accurate simulations for optimizing existing products or developing new ones, thereby shortening the production cycle.Figure 16 shows the predictive ability of different models for flow stress.The coefficient of determination  for the Arrhenius model is 0.977; for the HS model, it is 0.967; and for the new model, it is 0.99.Combined with the previous analysis results, the new model has the best fitting effect under various strain rate and temperature conditions, accurately predicting the stress-strain relationship of the material and having the highest predictive accuracy.The Arrhenius model follows, with good predictive ability, but it is slightly inferior to the new model.The HS model has the poorest predictive ability, with significant deviations, especially under high-stress conditions.Therefore, the new model demonstrates the highest reliability and accuracy in practical applications.Consequently, this high-precision model, with its mathematical expression, can be integrated into forming simulation software.This integration provides accurate simulations for optimizing existing products or developing new ones, thereby shortening the production cycle.

Conclusions
In this paper, a new constitutive model based on Taylor series expansions and partial derivatives is proposed.Hot compression tests were conducted at varying strain rates and temperatures to provide the necessary experimental data for model development.The new model significantly improves prediction accuracy without increasing the complexity of material parameters.The predictability of the new model and classical models (Arrhenius and HS models) was comprehensively analyzed, with mathematical proofs provided for the limitations of the classical models.The main conclusions are as follows: (1) The representative strain levels (low strain ε = 0.136, medium strain ε = 0.518, and high strain ε = 0.9) were selected.By analyzing the partial derivatives of logarithmic stress with respect to temperature and the different orders of logarithmic stress with respect to logarithmic strain rate, it was found that the relationship between logarithmic stress and temperature can be accurately described by a quadratic polynomial function, and the

Conclusions
In this paper, a new constitutive model based on Taylor series expansions and partial derivatives is proposed.Hot compression tests were conducted at varying strain rates and temperatures to provide the necessary experimental data for model development.The new model significantly improves prediction accuracy without increasing the complexity of material parameters.The predictability of the new model and classical models (Arrhenius and HS models) was comprehensively analyzed, with mathematical proofs provided for the limitations of the classical models.The main conclusions are as follows: (1) The representative strain levels (low strain ε = 0.136, medium strain ε = 0.518, and high strain ε = 0.9) were selected.By analyzing the partial derivatives of logarithmic stress with respect to temperature and the different orders of logarithmic stress with respect to logarithmic strain rate, it was found that the relationship between logarithmic stress
, respectively, show the comparison of the predicted and experimental flow stress values for the Arrhenius model and the HS model.It can be seen that the error is larger at lower temperatures of 850 • C, 900 • C, and 950 • C. As shown in Figure 6, the predictive accuracy of the HS constitutive equation is lower at lower temperatures (850 • C, 900 • C, and 950 • C), while it is higher in the high-temperature range (1000 • C, 1050 • C, 1100 • C, 1150 • C, and 1200 • C).

Figure 7 .
Figure 7. Different orders of partial derivatives of logarithmic stress with respect to temperature at a low strain level (0.136): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 7 .
Figure 7. Different orders of partial derivatives of logarithmic stress with respect to temperature at a low strain level (0.136): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 7 .
Figure 7. Different orders of partial derivatives of logarithmic stress with respect to temperature at a low strain level (0.136): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 8 .
Figure 8. Different orders of partial derivatives of logarithmic stress with respect to temperature at a medium strain level (0.518): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 8 . 25 Figure 9 .
Figure 8. Different orders of partial derivatives of logarithmic stress with respect to temperature at a medium strain level (0.518): (a) zeroth order; (b) first order; (c) second order; and (d) third order., x FOR PEER REVIEW 14 of 25

Figure 9 .
Figure 9. Different orders of partial derivatives of logarithmic stress with respect to temperature at a high strain level (0.9): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 10 .
Figure 10.Different orders of partial derivatives of logarithmic stress with respect to logarithmic strain rate at a low strain level (0.136): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 10 .
Figure 10.Different orders of partial derivatives of logarithmic stress with respect to logarithmic strain rate at a low strain level (0.136): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 11 .
Figure 11.Different orders of partial derivatives of logarithmic stress with respect to logarithmic strain rate at a medium strain level (0.518): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 12 .
Figure 12.Different orders of partial derivatives of logarithmic stress with respect to logarithmic strain rate at a high strain level (0.9): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 11 . 25 Figure 11 .
Figure 11.Different orders of partial derivatives of logarithmic stress with respect to logarithmic strain rate at a medium strain level (0.518): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 12 .
Figure 12.Different orders of partial derivatives of logarithmic stress with respect to logarithmic strain rate at a high strain level (0.9): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

Figure 12 .
Figure 12.Different orders of partial derivatives of logarithmic stress with respect to logarithmic strain rate at a high strain level (0.9): (a) zeroth order; (b) first order; (c) second order; and (d) third order.

k 10 ,
k 11 , and k 12 are material parameters.The constant term k 1 represents the reference value of stress with respect to the strain rate and temperature.The first-order term k 2 ln .

Materials 2024 ,
17,  x FOR PEER REVIEW 18 of 25 partial derivatives must be determined based on the complexity of their mechanical behavior.The accuracy of the new constitutive model is also ensured by these factors.Additionally, the coefficients must be solved mathematically.

Materials 2024 , 25 (
17, x FOR PEER REVIEW 21 of ) is 0.9948; the root mean square error (RMSE) is 22.5; the sum of squared errors (SSE) is 16,356; and the sum of absolute errors (SAE) is 5561 MPa.These values indicate that the model has superior fitting and predictive capabilities.The Arrhenius model follows closely, with a correlation coefficient () of 0.9888, a root mean square error (RMSE) of 44.1020, a sum of squared errors (SSE) of 19,450, and a sum of absolute errors (SAE) of 8832 MPa.However, its overall performance is slightly inferior to the new model.The HS model shows the poorest performance in all indicators, with a correlation coefficient () of 0.9827, a root mean square error (RMSE) of 44.0915, a sum of squared errors (SSE) of 38,881, and a sum of absolute errors (SAE) of 10,877 MPa, indicating its relatively weaker predictive ability.

, 17 ,
x FOR PEER REVIEW 22 of 25

Table 1 .
The normal chemical composition of the GH4169 superalloy (wt%).

Table 1 .
The normal chemical composition of the GH4169 superalloy (wt%).
Figure 1.Original microstructure of the wrought GH4169 superalloy.

Table 2 .
High-temperature isothermal compression test parameters.

Table 2 .
High-temperature isothermal compression test parameters.

Table 3 .
Stress matrix corresponding to different strain levels (the units of strain rate, temperature, and stress are s −1 , K, and MPa, respectively).

Table 4 .
Multivariate linear regression results of material parameters.

Table 5 .
Relationship between the coefficients of the new constitutive model for GH4169 alloy and strain.

Table 5 .
Relationship between the coefficients of the new constitutive model for GH4169 alloy and strain.
σi and σ i represent the predicted and experimental values of the ith data point, respectively, σ and σ represent the mean predicted and experimental values, respectively.The value of R ranges from [−1, 1].If R is closer to 1, it indicates a stronger positive correlation between the predicted and experimental values.If R is closer to −1, it indicates a stronger negative correlation, while values close to 0 indicate no significant linear correlation between them.Both the root mean square error (RMSE) and the sum of squared errors (SSE) measure the overall magnitude of the deviation between the model predictions and the experimental values.However, RMSE focuses more on reflecting local deviations, while SSE emphasizes reflecting global deviations.Their calculation formulas are shown in Formulas (

Table 6 .
Quantitative comparison of predictive accuracy of different models.